Optimizing the spin sensitivity of grain boundary junction nanoSQUIDs — towards 
detection of small spin systems with single-spin resolution 
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We present an optimization study for improving the spin sensitivity of nanoSQUIDs based on 
resistively sliunted grain boundary Josephson junctions. In addition the dc SQUIDs contain a 
narrow constriction onto which a small magnetic particle can be placed with efficient coupling of 
its stray magnetic field to the SQUID loop. The separation of the location of optimum coupling 
from the junctions allows for an independent optimization of coupling and junction properties. Our 
analysis of the dependence of the flux noise of the SQUID and of the coupling of a magnetic particle 
on the geometric parameters of the SQUID layout shows that a spin sensitivity of a few fiB /"v/Hz 
should be feasible for optimized parameters, respecting technological constraints. 

PACS numbers: 85.25.CP, 85.25.Dq, 74.78. Na, 74.72.-h 74.25.F- 74.40.De 
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I. INTRODUCTION 

Miniaturized direct current (dc) superconducting 
quantum interference devices (SQUIDs) with dimensions 
in the sub-micrometer range (nanoSQUIDs) arc promis- 
ing devices for the sensitive detection and investigation 
of small spin systems.— The basic idea behind this is to 
attach a small (nanometer-sized) magnetic particle di- 
rectly to the SQUID and trace out magnetic hysteresis 
loops of the particle. This shall be done by detecting the 
change of the stray magnetic field of the particle with 
magnetic moment fl via the change of the magnetic flux 
$ coupled to the SQUID loop<^"— . In order to meet the 
ultimate goal of detecting the flipping of only a few elec- 



tron spinS)^ the spin sensitivity S, 



1/2 



has to be 



optimized carefully via reducing the spectral density of 
flux noise 5$ of the SQUID and increasing the coupling 
factor 0^ = (with /i = |^|). can be reduced 

by shrinking the size of the SQUID loop, and hence its 
inductance L, and 0^ can be increased by placing the 
particle on a narrow constriction inserted in the SQUID 
loop, which motivates the need to implement sub-micron 
SQUID structures. 

Until now, the most common approach for the realiza- 
tion of nanoSQUIDs is to use constriction type Josephson 
junctions intersecting small SQUID loops (sec e.g. Ref. @ 
published in a special issue on nanoSQUIDs and related 
articles therein). However, this approach comes with 
several drawbacks: Constriction type Josephson junc- 
tions often show hysteretic current-voltage characteris- 
tics (IVCs). This does not allow to operate the SQUID 
magnetometers continuously, which is required for the 
investigation of the magnetization dynamics of the sam- 
ple under investigation. Furthermore, this does not allow 
to use advanced SQUID read-out electronics, which have 
been developed for SQUIDs with non-hysteretic Joseph- 
son junctions and which are required for the read-out of 
ultra-sensitive dc SQUIDs. Also, the noise properties of 
such junctions are not well understood and hence hard to 
optimize. And, finally, the magnetic particles have to be 



placed close to the constriction type junctions in order 
to achieve optimum coupling. However, this means that 
the junction properties and the coupling factor (p^^ can- 
not be optimized independently, which hampers a careful 
optimization of the spin sensitivity. 

Wc should note here, that very sensitive Nb thin film 
nanoSQUIDs based on constriction type junctions, resis- 
tively shunted with a thin W layer, have been realized.— 
However, in this case, the devices show optimum per- 
formance only in a narrow range of temperature T not 
too far below the transition temperature Tc of Nb, which 
makes them less interesting for applications. With re- 
spect to the application of nanoSQUIDs for the detection 
of the magnetization reversal of nanomagnets, the most 
interesting regime of operation is at T w 1 K and below 
and at very high magnetic fields in the tesla range ii It 
has been demonstrated that Nb thin film nanoSQUIDs 
based on constriction type junctions can be operated in 
impressive background fields up to 7T<^ However, the 
upper critical field of typical Nb thin films 1 T) re- 
quires to use very thin Nb films with thicknesses of only 
a few nm, i.e. well below the London penetration depth 
Xl of the Nb films, if such SQUIDs shall be operated in 
tesla fields. This leads to a large kinetic inductance con- 
tribution to the SQUID inductance, and hence a large 
fiux noise of such SQUIDs, which does not allow to use 
the huge potential for the realization of ultralow-noise 
nanoSQUIDs. 

In order to circumvent the above mentioned draw- 
backs, we recently started to develop dc nanoSQUIDs 
based on YBa2Cu307 (YBCO) thin films with submi- 
cron wide bicrystal grain boundary Josephson junctions 
(GBJs)<^ Due to the huge upper critical field of YBCO, 
such SQUIDs can be realized with film thicknesses on 
the order of Al and above and operated in tesla fields. 
Furthermore, due to the large critical current densities of 
the YBCO GBJs (several lO^A/cm^ at T = 4.2 K and 
below for a grain boundary misorientation angle of 24°) 
submicron junctions still yield reasonably large values of 
the critical current Iq . In order to achieve non-hysteretic 
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IVCs, the GBJs are shunted by a thin Au film. Due 
to the fact that the barrier of the GBJs is oriented per- 
pendicular to the YBCO thin film plane, it is possible 
to apply tesla magnetic fields in the plane of the film, 
without a significant reduction of /q'^ And finally, im- 
plementing an additional narrow constriction (which can 
be much narrower than the GBJs) in the SQUID loop, 
the optimization of the coupling factor is possible with- 
out affecting the junction properties. 

Here, we present a detailed optimization study of 
the spin sensitivity of such grain boundary junction 
nanoSQUIDs by analyzing the dependence of the flux 
noise 5$ and the coupling factor on the geometry of 
our devices. We find that for an optimized SQUID geom- 
etry a continuous detection of magnetic moments down 
to a spin sensitivity 5*^^ of a few ^.b/ \/Hz is feasible (/is 
is the Bohr magneton). 

II. NANOSQUID DESIGN 

The layout of the nanoSQUID (top view) is shown in 
Fig. [H The SQUID structure is patterned in a YBCO 
thin film of thickness d, covered by a thin Au film with 
thickness ^au- The two bridges straddling the grain 
boundary have a width wj and length Ij. The upper part 
of the SQUID loop contains a constriction of width Wc 
and length l^- An applied bias current / is flowing from 
top to bottom across the two GBJs. A small magnetic 
particle can be placed on top of the constriction, and 
an in-plane magnetic field (perpendicular to the grain 
boundary) can be applied without significant suppression 
of the critical current /q of the two GBJs. 

Optimizing the SQUID for spin sensitivity means to 
minimize the ratio S^/(f>j^. The coupling factor (/)^ is es- 
sentially determined by the geometry of the constriction, 
i.e., its width Wc and thickness d. 5*$ depends on the 
SQUID inductance L and on the junction parameters /q. 
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FIG. 1: (Color online) Schematic view of the nanoSQUID 
layout, divided (by white dotted lines) into the constriction 
(inductance Lc, length Ic, width Wc), two corners (each with 
inductance Le), the two junctions (each with inductance Lj, 
length Ij, width uij) and the bottom part (inductance !/(,)• 



resistance R and capacitance C. If the constriction could 
be made not only arbitrarily thin and narrow, but also 
arbitrarily short, one could envision a scenario, where 0^ 
reaches a value around O.S^oZ/iB^ while, at the same 
time, the inductance of the constriction remains small 
($0 is the magnetic fiux quantum). Then, 5*$ could be 
optimized independently by proper choice of the SQUID 
size and the junction properties. For the type of device 
we discuss here, this is certainly not the case and we thus 
look for an optimization, which is compatible with tech- 
nological limitations. A large coupling 0^ demands as 
narrow and thin as possible constrictions. On the other 
hand, for too narrow constrictions, given a fixed value of 
d, its inductance and thus also the total inductance L 
of the SQUID may become too large, possibly degrading 
the flux noise. This may be counterbalanced by choosing 
a different film thickness and changing, e.g., the junction 
width wj. 

In the following sections, we derive explicit expressions 
for the dependence of (Sec. IIIip and 5$ (Sec. lIVp on 
various geometric and electric SQUID parameters, which 
then allows us to optimize (Sec.|V|. 

III. COUPLING FACTOR 

Wc numerically calculate the coupling factor 0^ ~ 
<f>//i, i.e. the flux $ coupled into the SQUID loop by 
a point-like particle with magnetic moment /i, using 
the software package 3D-MLSI, which solves the London 
equationsiii Details on the calculation procedure can be 
found in Ref. 0. In brief, for a given SQUID geometry 
(cf. Fig. [T]), one calculates the magnetic field distribu- 
tion B{f) generated by a current J circulating aroimd 
the SQUID hole. The coupling factor is obtained from 
0p = —e^B{r)/J. Here, is the direction of the mag- 
netic moment /i at position r. 

For a particle with its magnetic moment oriented in 
the plane of the thin film SQUID perpendicular to the 
grain boundary, we find a maximum 0^ if the particle 
is placed directly on top of and as close as possible to 
the constriction.— Assuming that a particle is placed 
at a distance of 10 nm above the constriction (without 
an Au layer, which can be removed without affecting 
the junction properties), wc calculate 0^ in the range 
40 nm < Wc < 200 nm and 40 nm < d < 500 nm and find 
the dependence 

with do = 320 nm, wq = I nm and 7^ = 23n<I>o//iB. 

As expected, we find a decreasing (/)^ with increas- 
ing width Wc of the constriction. Our simulations yield 
0^ oc w""'^, i.e., a change in 0^ by a factor of 3 within the 
simulation range used for Wc- Furthermore, 0^ also de- 
creases with increasing film thickness d. However, within 
the (larger) simulation range for d, we find a much weaker 
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dependence on d; in particular, for d^dg, the coupling 
starts to saturate. We note that some uncertainty in de- 
termining (j)fj_ arises from the fact, that 3D-MLSI uses 
supercurrents that flow only within a 2D layer. Hence, 
in the constriction region - in particular in the case of a 
large aspect ratio d/wc ~ there may be a nontrivial 3D 
current flow. As a consequence, the prediction of a satu- 
ration of with increasing d may not be correct. On the 
other hand, we will show below, that 5$ shows clearly a 
stronger dependence on d than (/)^ does. Hence, for the 
optimization of 5*^ the uncertainty in is not very 

critical. 



IV. FLUX NOISE 

In order to determine the flux noise of the SQUID, we 
use the theoretical expression obtained from Langevin 
simulations, 



5$ = f{l3L)^okBTL/IoR, 



(2) 



which is valid for a Stewart-McCumber parameter 
= 27r/oi?2c/$o 5^ 1 and T^l < 0.1 Here, T = 
2'KkBT / Iq^q is the noise parameter, and /3l = 2i/o/$o 
is the screening parameter. For > 0.4, f{PL) ~ 
4(1 -I- /3l). For lower values of increases. 

The first factor to be discussed is IqR. The junction 
resistance R can be varied to some extent by varying the 
thickness Jau of the Au layer covering the YBCO film; 
the maximum achievable value is the unshuntcd junc- 
tion normal state resistance Rn (for ^au = 0). For 24° 
YBCO grain boundary junctions, IqRn values ~ 2— 3mV 
are achievable at 4.2 K.— However, such junctions typi- 
cally have hysteretic IVCs. We thus demand /3c ~ 1 to 
avoid hysteresis. Ideally, one would like to derive an 
expression for IqR as a function of wj, d and c^au us- 
ing the constraint /3c ~1 and assuming certain values for 
the critical current density joi unshunted normal junc- 
tion resistance times area p = R^wjd and capacitance 
per junction area C". However, the scaling of R with 
wj, d and c^au is currently not known. Furthermore, an 
estimate of C as a function of wj and d, based on vari- 
ous scaling laws available in the literature,—""—, is quite 
difficult, in particular since it is difficult to determine C 
for underdamped YBCO GBJs and since the stray ca- 
pacitance due to the commonly used SrTiOa substrates 
may play an important roleJ^ On the other hand, we 
have fabricated 24° YBCO GBJs with different junction 
width and film thickness {wj — 450 nm, d = 50 nm and 
wj = 330 nm, d ~ 300 nm), using the focused ion beam 
(FIB) milling technique as described ini^. Those junc- 
tions had only weakly hysteretic IVCs but IqR k, 1 mV 
at T = 4.2 K. Below we will find an optimum junction 
width wj ~ 260 nm (for d = 300 nm), which is not much 
below the width of those junctions. Thus, rather than 
introducing an ill-defined scaling of /q-R with wj and d, 
below we fix I^R and use 1 mV as a realistic value. 



We next determine the dependence of the SQUID in- 
ductance L on the various geometrical parameters. We 
separate the SQUID into the constriction (inductance Lc, 
length Zc, width Wc), the two (symmetric) bridges con- 
taining the junctions (inductance ij, length Ij, width 
wj), the two corners connecting the constriction and the 
junction arms (inductance ig), and the bottom part of 
the SQUID (inductance L;,), as indicated in Fig.[T] Then, 
L is given by 



L = Lc + 2Lj + 2Le 



(3) 



We should find LdwcJcd), Lj{wj,lj,d), Le{wc,wj,d) 
and Lb{lcTWj,d). Trivially, Ic and Ij should be as small 
as possible. From simulations, using 3D-MLSI with = 
335 nm^ we find the parametrization Lc{wc,lc,d) « 
L' ■ lc/{wcd), where L' = 143pH-nm. This expression 
fits the simulated Lc well, within the parameter range 
40 nm < Wc < 300 nm and 40 nm < d < 500 nm. cov- 
ered by the simulations. We use the same parametriza- 
tion for Lj(wj,lj,d). For the corners we find, within a 
15% variation with respect to wj and Wc, the expres- 
sion Le « I^'e/d, with Lg = lOOpH-nm. Finally, we 
find Lb « L'fjlc/wjd + L'^/d, with = 45pH-nm and 
L'll = 150pH-nm. Inserting these expressions into Eq. ([3]) 
yields 



L ; 



~d 



{2lj + Uc) 



wj 



(4) 



with b = L'JL' = 0.315 and r = (2L^ -I- L'l)/L' = 2.45. 

For the minimization of S*^, we will use /3l as a variable 
parameter. Since both, L and wj arc not independent 
of each other and arc related to /3l, we express both as 
functions of 13 l. This will allow us to eliminate L and 
wj in the final expression for which has to be opti- 
mized. With /3l = 2/oL/$o and Iq = jowjd, we obtain 
wj{(3l,L) = <I>o/^L/2jodi. Inserting this into Eq. (g]) 
yields 

(5) 



d \Wn 



K 



with k{Ij, IcJo) = '2(21 J + blc)jaL'/^o- 

Inserting Eq. © into Eq. ^ and using /(/3i) = 4(1 - 
Pl) finally yields 



c ^ e do f lc 
d \Wc 



1 + 

^ I3l 



(6) 



with S^l = 2^1^^^ = 7.1n$o/%/H^ at T = 4.2 K 
and with IqR = 1 mV. The most important result here 
is the scaling S** cx 1/d. This is due to the fact that the 
SQUID inductance L cc 1/d within the simulation range 
for d because of the increase of the kinetic inductance 
contribution with decreasing d below Al. For dZ2\L we 
expect a saturation of L{d) and hence of S^{d). However, 
we will neglect this for the optimization of S*^, since val- 
ues for dz500 nm are outside the simulation range for 
L(d) and since we cannot expected to produce high- 
quality GBJs for such large values of d. 
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FIG. 2: (Color online) Spin sensitivity S]J'^{I3l), calculated 
from Eq. (O at T = 4.2 K for 10 values of film thickness d 
(in 50 nm steps), with constriction width ~ 50 nm (and 
Ic = 200 nm, /,/ = 400 nm, jo = 3mA//im^). 



V. OPTIMIZATION OF SPIN SENSITIVITY 
VIA IMPROVED SQUID GEOMETRY 

With Eq. (m and ([1]) wc can express the spin sensitivity 

as 

al/2 ^ cl/2 "/i ( lie j l + h 

^ -^,,0^2.3^ + 1 Uoj V^c+ V^~^ 

(7) 

with Sl(^ EE 5^^0/7^ = 0.3lMB/\/Hi at T = 4.2K and 
with IqR = I mV. 

The evaluation of Eq. ([7]) shows that the minimiza- 

1/2 

tion of Sfj. {f3L,d,Wc,lc,lj) requires to make Ic and Ij 
as small as possible. Due to technological limitations, we 
use Ic = 200 nm and Ij = 400 nm below. The constriction 
width enters Eq. ([7]) via the factors w'^''' {Ic/wc + J")^^^, 

demanding Wc to be made as small as possible. Fur- 

1/2 

thermore, decreases with increasing d within the 

simulation range; for large enough values of d, we find 
cx l/^fd. Hence, the film thickness should be as 
large as possible, respecting technological constraints. As 
mentioned at the end of Sec. IIVI we expect L and hence 
5$ to saturate for d i 2Al. In this case, which lies outside 
our simulation range, we expect to increase again 
with increasing d, due to the degradation of (f)^. Hence, 

1/2 

we should expect a minimum in S^, for an optimum d, 

which however lies above our simulation range and also 

above technologically feasible values for d. 
1/2 

The dependence 5*^ {13 l) at T = 4.2K is shown in 
Fig- m for various values of d, with fixed Wc = 50 nm. All 
curves exhibit a shallow minimum at a value which can 
be calculated from Eq. (O as /Ji^miTi = ^(1 + Vl + k^^). 
With the (fixed) values for Ic and Ij and with jo = 
3mA//xm^, we obtain k = 0.36 and PL.min = 1.05. For 
a reasonable value d — 300 nm (cf. thick line in Fig. [2]) 



1/2 / — 
we obtain ~ 33n$o/vHz, (j)^ = 20n<l5o//^B and 

~ 1.6^b/a/Hz. Corresponding SQUID parameters 
are listed in Table H] ("opt. device 1"). 

This is a drastic improvement, compared to the ex- 
perimental value S]/"^ « GO^b/VHz for our, so far, 
best experimental device with parameters 0? = 60 nm, 
Wc = 90 nm, wj = 130 nm, /3i = 0.65, L ~ 36 pH, I^R = 

0.13 mV, = 1.3Ai$o/Vlfe and 0^ = 21n$o/MB'^ 

1/2 1/2 
The main improvement of Sp! comes from . Before 

discussing the various factors improving this quantity, 

we should say, however, that standard YBCO SQUIDs 

1 /2 

often have a , which is a factor of 2-10 higher than 

the theoretical predictionai^ (for the experimental device 

1/2 

mentioned above was a factor ~ 6 higher than pre- 
dicted by Eq. ([2])). In this respect we expect the theo- 

1 /2 

retical to be by a factor ^ 2-10 too low, even if all 
factors entering S*^^^ are feasible. 

1/2 

One major improvement of arises from the to- 
tal SQUID inductance L, which is almost an order of 
magnitude lower than the inductance of our experimen- 
tal device. This is mainly because of the much larger 
value of d. The low value of L seems realistic; in terms of 
geometrical factors, Wc = 50 nm and lengths Ic and Ij are 
chosen sportive, but not impossible. This is supported by 
preliminary tests, which show that 300 nm thick YBCO 
films can be patterned properly by FIB milling. 

1 /2 

Another big improvement of comes from the 

optimized value of IqR^ which is by a factor of ^ 8 
higher than that of our experimental device. For nar- 
row (wj ~ 100 nm) and thin {d ^ 50 nm) junctions 
this would be impossible to achievej2. However, the opti- 
mization yielded much wider (wj 260 nm) and thicker 
(d ~ 300 nm) junctions, for which we have already found 
that IqR ^ 1 mV is achievable. 

Finally, we note that, if we take more easily achiev- 
able values d = 200 nm and Wc = 100 nm (other input 
parameters are the same as for the initial optimization), 
we still get 5"^^ = 2.6^b/-\/Hz, with parameter values 
wj = 375 nm, L = 4.8 pH, Iq = 225 //A and R = 4.44 
(see Table. H] for parameters of "opt. device 2"). The re- 

1 /2 

suiting is less than a factor of 2 different from the 
value obtained for the initial input parameters, so that 
we are optimistic that the design discussed in this paper 
is indeed capable of obtaining a spin sensitivity in the 
range of a few /is/VHz. 

VI. CONCLUSIONS 

In summary, we have performed a detailed analysis 
of the coupling factor and the spectral density of 

1 /2 

flux noise S'$, and hence of the spin sensitivity S^, = 

1 /2 

/(f>^ of a grain boundary junction dc nanoSQUID. 
Based on the calculation of (j)^ and 5*$ , we derived an ex- 
plicit expression for the spin sensitivity 5"^^ as a function 
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units 


d 
nm 


nm 


nm 


nm 


Wj 

nm 


Pl 


L lo 
pH fiA 


R 

n 


loR 
mV 


jo Lc Lj 
mA//im^ pH pH 


pH 


Li 
pH 


ol/2 

n$o /"v/Hz 


(t>tJ. 
n<l>o / fiB 


q1/2 
flB 1 \/Hz 


opt. 
device 1 


300 


200 


400 


50 


260 


1.05 


4.7 233 


4.30 


1 


3 1.9 0.74 


0.33 


0.62 


32.7 


19.9 


1.65 


opt. 

device2 


200 


200 


400 


100 


375 


1.05 


4.8 225 


4.44 


1 


3 1.4 0.76 


0.50 


0.87 


33.3 


12.7 


2.62 



TABLE I: Summary of geometric and electric nanoSQUID parameters (as defined in the text). The values for "opt. device 
1" are calculated for an optimized device with film thickness d — 300 nm and constriction width Wc — 50 nm. The values for 
"opt. device 2" are calculated for more relaxed values d = 200 nm and Wc = 100 nm and otherwise identical input parameters 
for lc, Ij, IqR and (optimum) /Jl. 



of the geometric and electrical parameters of our devices. 

1/2 

This allows for an optimization of 5^ , which predicts a 
spin sensitivity of a few /is/VHz. Such a low value for 

1/2 

S)j, can be achieved by realization of very low induc- 
tance nanoSQUIDs with ultra-low flux noise on the order 
of only a few tens of n<I>o / V^Hz. This poses severe chal- 
lenges on proper readout electronics for such SQUIDs. 
It remains to be shown whether or not the readout of 
such ultralow-noise SQUIDs is feasible and whether or 
not the envisaged values for the spin sensitivity can also 
be achieved in high fields, which is a major driving force 



for using these grain boundary junction nanoSQUIDs. 
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